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ABSTRACT 

Hitherto, six P-type planets are found around five binary systems, i.e. Kepler-16 b, 34 b, 35 b, 38 
b, 47 b, c, which are all Neptune or Jupiter-like planets. The stability of planets and the habitable 
zones are influenced by the gravitational and radiative perturbations of binary companions. In this 
Letter, we check the stability of an additional habitable Earth-mass planet in each system. Based on 
our simulations in 10 Myr, a habitable "Earth" is hardly stable in Kepler-16 while a stable "Earth" 
in Kepler-47 close to the boundaries of the habitable zone is possible. In contrast, Kepler-34, 35 and 
38 seem to have high probabilities of being able to tolerate a stable "Earth" in their habitable zones. 
The affects of transit time variations are quite small due to the small mass of an undetected " Earth" , 
except that of Kcpler-16 b. With a time precision of 10"'^ day(^ 88s), an "Earth" in the corotational 
resonance with Kepler-16 b can be detected in 3 years, while habitable "Earths" in Kepler-34 and 
38 systems can be detected in 10 years. Habitable "Earths" in Kepler-35 and 47 are not likely to be 
detected in 10 years under this precision. 
Subject headings: stars: individual: (Kepler-16, 34, 35, 38, 47) 



1. INTRODUCTION 

Searches for Earth-like exoplanets, especially the Ke- 
pler mission, make great contributions to our knowl- 
edge about terrestrial planets. To date, about 66 Earth- 
like (< IOM0) planets are confirmecQ and more than 1100 
candidates with radii i?p < 2R^ are detected^ Con- 
sidering these Earth-like planets with minimum masses 
< lOAf0, only few of them have a moderate equilibrium 
temperature and are considered as habi table planets 
(IKopparapu et al..ir2013). e.g.. HD40307 g (iTuomi et al..l 



I2012D and GJ667C c (|Delfosse et al. .112013) . Meanwhile, 
11 candidates with radii < 2Rq are in the habitable zone 
()Gaidos.ll2013D . Habitable Earth-like planets are likely 
to stay in more compacted multi-planet systems. For ex- 
ample, HD40307 contains six planets and GJ667C has at 
least two planets. 

Planets in habitable zone (hereafter HZ) must have 
a moderate temperature to maintain liquid water on 
the surfa ce arid hav e a possibilit y of sustaining life 
(JKasting eTaTI 119931: IJones et all '12001'). The tem- 
perature of a planet is usu ally considered via sev- 
eral conditions (jGaidos .1120131) . e.g ., planetary albedo, 
orbit al eccentricity and obliquity (jWilliams &: Pollard.1 
I2003D. components of the planetary atmosphere 



( Pierrehumbert et al. .1120111 ) . etc. In single star systems, 
the HZ is fixed due to the stable flux of the host star 
in main sequences. However, the HZ in binary systems 
is different due t o the radiative pe rturbation of the bi- 
nary companion. lEggl et all ()2012[ ) provided an analytic 
method to estimate the HZs for S-type planetary orbits, 
and pointed out that HZs are strongly affected by the 
eccentricity of the binary. Kane & Hinkel (2013) stud- 
ied HZs for P-type planetary orbits, which would oscil- 
late with time due to the different orbital phases of the 



binary stars. The stabilities of planets in binary sys- 
tems are also different fro m those in single star systems 
(jHolman &: Wieger^ 119991 ). Considering a multi-planet 
system, planets with equal masses are stable in single 
star systems if their separations exceed a c ritical value 
~7 - 10 HiU radii (jChambers et al..lll996t iZhou et al..l 
l2007f l. Adding the perturbation of a binary companion, 
the stabilities of planets around the binary are queried. 
In this Letter, we foc us on the five circum- 
binary systems: Ke pler-1 6 (iDovle et al.,|[20Tll) . 34, 35 
([Welsh et al.J [2011) 38 (lOroze et al..l l2012bD and 47 
(lOroze et al..l l2012aD . All six planets detected are 



Jupiter-like or Neptune-like planets(see Table 1), which 
probably have a thick atmosphere and may not form a 
sol id surface, thus the y are not suitable for life. As noted 
bv lGong et al..l (|2013[ ). the formation of habitable Earth- 
like planets in circumbinary systems is possible. Here 
we assume that an additional "Earth" formed in these 
systems to obtain the tolerant regions of stable habit- 
able "Earths". We also simulate the transit time varia- 
tions (hereafter TTVs) of the existing planets with and 
without an additional "Earth". The differences between 
these two TTVs may help us to predict the existence 
of an undetected "Earth" through fur ther observations 
(|Agol et al.ll2005l: iLithwick et al. .112011) . 

2. ADDITION-" EARTH" MODEL 

Researchers usually estimate the boundaries of HZ via 
the maximum greenhouse effects provided by a CO2 at- 
mosphere and via runa way greenhouse effect (i.e., lo ss 
of water. (iKasting et a ll993t lUnd erwood et alJ 2003|)). 



As IKopparapu et al.J (|2013D pointed out most recently, 
the stellar fluxes (with unit 1368 W m~^) at the bound- 
aries of HZ are 



^ Data from 'http://exoplanets. orgl 
^ Data from http://www.nasa.govl 
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where T* — Teff — 5780K and the va l ue oi S^s^p^^a^h^c 
and d are Hsted in iKopparapu et all (|2013D . see Tab. 3 
therein. In single star systems, Tos is the effective tem- 
perature of the star. In binary systems, the equivalent 
Toff of the binary are calculated by Wien's displacement 
law after combining the stell ar spectral radiances o f bi- 
nary stars (see Section 3.2 in iKane fc Hinkell ()2013[ )). 

The equilibrium temperature of a planet Tp is de- 
termined in two steps. (l)Estimate the combined flux 
on the planet S \iy S = S\ -\- 5*2, where Si = (t(1 — 
a){Ri/ri)'^T^^ ;, (i=l,2), where a is the Boltzmann con- 
stant, a is the albedo of the planets, Ri is the radius of 
star J, r^ is the distance between the planet and star i and 
Teff_i is the effective temperature of star i. (2)Estimate Tp 
of the planet by energy equilibrium: Trr^ x S — Airar'iT^, 
therefore we obtain 



Tp = {S/Aa) 



1/4 



(3) 



Substituting Equations ([Ij and ^ into Equations ([3]) , we 
get the temperatures at the HZ boundary. Tin and Tout- 
Note that Tjn and Tout are independent of albedo a, but 
the location of the HZ depends on a in the expression of 
5*1 and 82- 

The configurations of the five observed circumbinary 
systems are listed in Table 1 according to observations. 
We simulate 100 cases for each system. In each case, 
we put an additional "Earth" with a mass= IM^ and 
let the "Earth" evolve under the gravities of both bi- 
naries and the existing planets. All these "Earths" are 
put in circular orbits with an inclination=l° initially; the 
phase angles of the " Earths" are chosen randomly. The 
initial locations of these "Earths" in each systems are 
interpreted in Section 3. Using the MERC URY package 
for TV-body simulations ([Chambers! 119990 . we estimate 
their stabilities and habitats based on our simulations. 
Due to the variation of Tp , only those "Earths" with 
Tp[t) G [Tout, Tin] are considered habitable in this work. 

3. RESULTS IN FIVE SYSTEMS 
3.1. Kepler- 16 

The equivalent temperature of binaries in the Kepler- 
16 system is Toff — 4307. 8K, and the HZ with tempera- 
tures [200.2K,271.6K] (Figure 1(a)) is located from 0.4 to 
0.75 AU approximately. The locations of HZ boundaries 
change little with time because of the small contribu- 
tion to the combined flux from the much cooler com- 
panion star. Note that the locations of HZ boundaries 
here are averaged by different azimuths. According to 
IQuarleset aUlftFm . "Earths" inside 0.66 AU (« 222K) 
are unstable in the system. The initial locations of our 
100 "Earths" are chosen between [173K, 222K] randomly. 

The red/blue circles represent the initial locations of 
unstable/stable "Earths" during our simulations (the 
same notations are used in all figures hereafter). Af- 
ter 10 Myr evolution, more than 90% of "Earths" are 
ejected mainly due to the perturbations of Kepler-16 b, 
and only seven "Earths" survive outside the HZ. The or- 
bits of other detected giant planets are also allowed to 
evolve under the gravities of planets and stars in our sim- 
ulations. However, due to the large mass of Kepler-16 b, 
if the "Earth" is stable, the orbits of the presently ob- 
served planets are nearly unchanged. The same is true 
for other systems in Section 3. Finally, we obtain a null 



"Earth" in the HZ of the Kepler-16 systems. As shown 
in Figure 1(b), we find two stable regions: r ^^0.72 and 
0.95AU. The inner two are in the corotational resonance 
with Kepler-16 b. Averaged by the time, their mean 
Tp w207.5 and 207.4 K, with variations about 19.7 and 
16.2 K. Therefore, these two "Earths" could be in the 
HZ only by considering some other different atmosphere 
models. 

We set the albedo a = in Figure 1. If we change the 
albedo to a=0.3 comparable with the Earth and a higher 
value 0.5 to model metal-rich planets with a smooth sur- 
face, the equilibrium temperature of planets Tp becomes 
0.915 and 0.841 times of that with a = 0. The two plan- 
ets in the same orbit with Kepler-16 b will obtain a lower 
Tp ~ 189. 8K for a = 0.3 and Tp - 174.5K for a = 0.5, 
which no longer seems likely to be habitable. 

3.2. Kepler-Si 

The Toff of Kepler-34 binary is 5892 K, and an "Earth" 
with Tp in [214.2K, 282.9K] is habitable (Figure 2(a)). 
The HZ is about 1.6-2.8 AU according to the fitting curve 
of Tp in Figure 2(b) (hereafter the same). The HZ bound- 
aries vary with time more obviously than others because 
of the large eccentricity of binary stars and their similar 
Toff. However, the variation of the inner HZ boundary 
is still very small(< 0.02 AU) since the HZ is far away 
from the binary compared with the separation of bina- 
ries. The initial locations of our 100 "Earths" are chosen 
between [173K, 570K] randomly. 

The green filled circles represent the initial locations of 
stable habitable " Earths" (hereafter the same) . Compar- 
ing with Kepler-16, the HZ in Kepler-34 is more stable 
for an additional "Earth". As shown in Figure 2(a), 74 
"Earths" far away from Kepler-34 b easily survive af- 
ter 10 Myr including 32 (75%) of the initial habitable 
"Earths" that remained in the HZ throughout that time. 
Figure 2(b) shows the locations of the 74 stable "Earths" 
and their Tp. The Tp{r) is similar to that in single star 

systems as a power law: Tp = 359.1K(r/lAU)~^/^. We 
find a stable boundary of about 1.5AU, outside which 
"Earths" are stable. The mean distance of the outer HZ 
boundary is about 2.8 AU. Thus, the Kepler-34 system is 
likely to tolerate a habitable "Earth" between 1.6 and 2.8 
AU when a = 0. Due to the perturbation of Kepler-34 
b, "Earths" close to it have a larger variations of Tp. 

Changing the albedo a to 0.3 and 0.5, the location of 
the HZ changed to an inner region, thus the outer bound- 
ary of HZ becomes 2.34 AU for a = 0.3 and 1.98 AU for 
a = 0.5. Overlapping with the stable boundary >1.5 
AU, we conclude that the tolerant region of habitable 
Earths in the Kepler-34 system becomes: 1.6-2.34 AU 
with a = 0.3 and 1.6-1.98 AU with a = 0.5(see Table 2). 

3.3. Kepler-35 

The equivalent temperature of binaries in Kepler-35 is 
Toff = 5435.5 K, and the HZ is from 1.2 to 2.0 AU with 
temperature [210.0K, 278.9K] (Figures 2(c) and (d)). 
Both the separation and eccentricity of the binary or- 
bit are smaller than that of Kepler-34, therefore, the HZ 
boundaries are less influenced(< 0.01 AU) by the mo- 
tion of binaries with time. The same effect occurs in 
Kepler-38 and 47, in which the HZ boundaries are nearly 
circular. The initial 100 "Earths" are choosen between 
[173K, 520K] randomly. 



Figure 2(c) indicates that 15 "Earths" in the in- 
ner region are unstable while other 85 "Earths" out- 
side still orbit the binary stars after 10 Myr, includ- 
ing 36 (-90%) of the initial habitable "Earths" that 
remained in the HZ. The Tp of the 85 stable "Earths" 
are shown in Figure 2(d), well fitted by a power law 
Tp = 298.8K(r/lAU)-i/2. We find that the tolerant re- 
gion of a stable habitable "Earth" with a = in the 
Kepler-35 system is about 1.2-2.0 AU and the variations 
of Tp of "Earths" in the HZ are all less than 8 K. 

Taking q:=0.3 and 0.5, the location of the HZ will move 
to an inner region, i.e., 1.0-1.67 AU for a = 0.3 and 0.85- 
1.41 AU for a — 0.5. Because smaller HZ overlaps with 
part of the unstable region, a smaller fraction of surviving 
"Earths" remain in the HZ in comparison with a — 0(see 
Table 2). 

3.4. Kepler-38 

With a radius of 0.3964i?j, a mass of - 0.078Mj for 
Kepler-38 b is estimated by the density of Neptune in 
our simulations. The equivalent temperature of binaries 
is Teff = 5496. 7K, and the HZ boundary is [210.5 K, 279.4 
K], about 1.6-2.9 AU (Figures 2(e) and (f)). Here The 
initial 100 "Earths" are chosen randomly from 1.2-4.3 
AU with temperature [173K, 350K]. 

Due to the tight binary orbit, the large distance be- 
tween "Earths" and the small Kepler-38 b, "Earths" in 
this system are all stable after 10 Myr (Fig.2e). Fi- 
nally, 38 (>90%) of the initial habitable "Earths" re- 
maining in the HZ. Fig.2f shows the fitting curve Tp — 

359.7K(r/lAU)^^/^ and indicates that the boundary of 
the HZ is from 1.6 to 2.9 AU with a = 0. The mean 
variations of Tp of "Earth" in the HZ are all less than 
lOK. 

The HZ becomes 1.34-2.43 AU for a = 0.3 and 1.13- 
2.05 AU for a = 0.5. Since all 100 "Earths" in the 
Kepler-38 system are stable, most "Earths" in the HZ 
initially are still habitable due to our simulations while 
changing the albedo to 0.3 or 0.5(see Table 2). 

3.5. Kepler- ^7 

Using the density of Neptune, we estimate the masses 
of Kepler-47 b, c as 0.0248Mj and 0.0933M,/, respec- 
tively. Their eccentricities are only limited as < 0.035 
and < 0.411, respectively. We choose a mean value for 
each planet, i.e., 0.035/2 ~ 0.017 for 47b and 0.411/2 - 
0.205 for 47c in our simulations. The equivalent tempera- 
ture of binaries is T^g = 5501. 5K, and the HZ boundary 
becomes [210.6 K, 279.4 K], about 0.9-1.6 AU (Figure 
2(g), (h)). The initial 100 "Earths" are chosen randomly 
between [173 K, 550 K], about 0.2-2.6 AU. 

After 10 Myr evolution, only 66 "Earths" survive and 
7 of them are habitable, as shown in Figure 2(g). Since 
Kepler-47 c is located in the HZ with eccentricity 0.2, it 
sweeps out a region of 0.8-1.2 AU, as shown in Figure 
2(h). Similar to the Kepler-16 system, Kepler-47 c also 
captures an "Earth" in the 1:1 orbital resonance with a 
large eccentricity ~ 0.65; therefore, this "Earth" has a 
quite large variation of Tp {^ 220K) and is not habit- 
able. Only "Earths" from 1.2 AU to 1.6 AU survive as 
habitable "Earths". Due to the perturbation of Kepler- 
47 c, these habitable "Earths" have a large dTp --20K. 
Inside the inner boundary of the HZ, there are 14 stable 



"Earths", with an even larger dTp > 30 K due to the 
perturbations of both Kepler-47 b and c. 

If the albedo a = 0.3, the HZ extends from 0.75 to 1.34 
AU, and there will be two tolerant regions for a stable 
habitable "Earth": 0.75-0.8 AU and 1.2-1.34 AU. With 
a — 0.5, the HZ becomes 0.64-1.13 AU, and the tolerant 
region becomes 0.64-0.8 AU(see Table 2). 

4. TTVS VIA THE ADDITIONAL "EARTH" 

Because all the "Earths" have an inclination=l°, we 
can't observe their transits. To determine the existence 
of stable "Earths" in the HZ, we can calculate the tran- 
sit time (Ttr) and TTVs of the existing planet in each 
system, which can be observed by follow-ups. We calcu- 
late two Ttr with and without an additional "Earth" in 
each systems, and check their differences dTt^- The ini- 
tial conditions used here are according to our results in 
Section 3. The "Earths" labeled as "sample" in Figure 1 
and 2, are chosen as the sample cases. The initial condi- 
tions of both binary stars and planets in this section are 
the same with the configurations in the sample cases at 
the end of simulations (time=10 Myr) in Section 3. 

Figure 3 shows the differences of both transit time 
(dTtr) and TTVs (dTTV) in 4000 days in each system. 
Note that we only show the transits of the larger bi- 
nary star A in each system. For Kepler-47 system, only 
the transits of 47 b are considered because the period of 
47 c is so long that few transit events can be observed. 
In Figure 3(a), all systems have increasing variations 
dTtr. In Kcpler-16, the "Earth" is in the corotational 
resonance with Kcplcr-16 b, and the interaction between 
these two planets is strong. Therefore, the dTtr is more 
obvious than other four systems. The dTtr of Kepler-16 
b can achieve ~ 10^^ day in 3 years and ~ 1 day in 10 
years. The other four systems have small variations due 
to the limited perturbation of the additional "Earth": 
dTtr < 10~^ day in 3 years and dTt^ < 10~^ day in 10 
years. We can sort these five systems by dTtr'. Kepler- 
16 b>Kepler-34 b~Kepler-38 b>Kepler-35 b>Kepler-47 
b. Figure 3(b) also shows the calculation of the TTV of 
the existing planets. Similar to dTtr, the strong interac- 
tion between Kepler-16 b and co-orbit "Earth" results in 
a larger dTTV'^ 0.1 day. The other four systems have 
smaU dTTV< 10^^ ^j^y. 

5. CONCLUSIONS 

With the motivation of finding additional "Earths" in 
the HZs of Kepler-16, 34, 35, 38 and 47, we simulated 
100 cases in each system and show the stability and hab- 
itability of these "Earths" (see Section 3). After 10 Myr 
evolution, most "Earths" are unstable in Kcplcr-16 and 
we do not find any habitable "Earth". However, two 
"Earths" in the corotational resonance with Kepler-16 b 
survived near the outer HZ boundary when the albedo 
a ~ 0. In Kepler-47, due to the perturbation of Kepler- 
47 c, only eight "Earths" from 1.2 to 1.6 AU, survive and 
are habitable. In Kepler-34, 35 and 38, "Earths" in the 
HZ are stable and most "Earths" in the HZ initially can 
remain in the HZ. The tolerant region of an additional 
habitable "Earth" in each system is shown in Table 2. 

Different albedos {a —0.3 and 0.5) of the "Earth" make 
the conclusions a little different. Comparing with a = 0, 
"Earths" absorb less flux and have a lower Tp, i.e., the HZ 



migrates to a closer and narrower region. In the Kepler- 
16 system, we still do not obtain any habitable "Earth" 
with a=0.3 or 0.5. Kepler-34, 35, and 38 systems have 
fewer habitable "Earths" that survived and the outer 
boundaries of the tolerant regions shrink. Overlapping 
the HZ with the unstable region in Kepler-47, there are 
two tolerant regions for a habitable " Earth" in Kepler-47 
when we set a = 0.3, as shown in Table 2. 

Comparing the Ttr with and without the perturbation 
of the additional "Earth", we show that: an "Earth" in 
the same orbit with Kepler-16 b can influence the Ttr of 



Kepler-16 b obviously, with dTt,- - 1 day and dTTV- 0.1 
day in 10 years. Meanwhile, dTtr and dTTV due to the 
perturbation of the habitable "Earth" in other systems 
are relatively small. With the precision of lO^'^ day ~ 
88s, we can identify dTTV of Kepler-16 b in three years. 
After 10 years of observations, the dTTV of Kepler-34 b 
and 38 b can be detected, while the dTTV of Kcplcr-35 
b and 47 b are too small to be detected. 

The work is supported by National Basic Research Pro- 
gram of China (2013CB834900) and National Natural 
Science Foundations of China (Nos. 10833001, 10925313, 
11003010). 
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TABLE 1 

Stellar and planetary characters in the five binary systems. 







Kepler-16 


Kepler-34 


Kepler-35 


Kepler-38 


Kepler-47^ 




Mass(Af0) 


0.6897 


1.0479 


0.8877 


0.949 


1.043 


Star A 


Radius{ijQ) 


0.6489 


1.1618 


1.0284 


1.757 


0.964 




Temperature(K) 


4450 


5913 


5606 


5623 


5636 




Mass(A/0) 


0.20255 


1.0208 


0.8094 


0.249 


0.362 


Star B 


Radius{ijQ) 


0.22623 


1.0927 


0.7861 


0.2724 


0.3506 




Temperature(K) 


2800 


5867 


5202 


3315 


3357 



semi-major axis{AU) 
binary period (days) 
parameters eccentricity 
rcfl-(Kli 



0.22431 

41.079220 

0.15944 

4307.8 



0.2288 

27.7958103 

0.52 

5892.0 



0.17617 

20.733667 

0.1421 

5435.5 



0.1469 

18.79537 

0.1032 

5496.7 



0.0836 

7.44837695 

0.0234 

5501.5 



Mass(Af,j) 0.333 0.22 0.127 0.3964ijj 0.27iJ,j 

semi-major axis(AU) 0.7048 1.0896 0.603 0.4644 0.2956 

Planet b period(days) 228.776 288.822 131.455 105.595 49.514 

eccentricity 0.0069 0.182 0.042 0.016 < 0.035 

Tp{a = 0)^ 200-215 310-390 370-420 460-490 480-510 



Mass(Mj) 
semi-major axis(AU) 
Planet c period(days) 

eccentricity 
Tp{a = 0) 



0.42i?j 
0.989 
303.158 
< 0.411 
267-273 



^Data are from lDovle et aI7]|2011l , F 



^Tgff is the effective temperature of the binary system 

'^Tp is the equilibrium temperature of planet with albedo a = 0. 



T1 [20Tg . lOroze et ani2012bl , lUroze et ani2012fal and lOroze et ai..|[20T2al respectively 



TABLE 2 

Main results in the five binary systems. 







Kcplcr-16 


Kcplcr-34 


Kepler-35 


Kcpler-38 


Kepler-47 


Stable "Earth" 




7 


74 


85 


100 


66 


Stable "Earth" 


in HZ 


<2 


32 


36 


38 


7 


Stable region (AU) 


>0.66 


>1.5 


>1 


>1.1 


0.4-0.8 , >1.2 


Tolerant region 


(AU) 












Q = 




~0.72 


1.6-2.8 


1.2-2.0 


1.6-2.9 


1.2-1.6 


a=0.3 




Null 


1.6-2.34 


1.0-1.67 


1.34-2.43 
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Fig. 1. — "Earths" in Kepler-16: (a) two green dotted lines are the HZ boundaries with albedo a = 0. Temperatures of each line are 
labeled. The two red dotted lines represent the extension of our 100 "Earths", with initial random locations shown in the plot. The 
red/blue circles represent unstable/stable "Earths" after 10 Myr evolution. The one labeled as "sample" is chosen for the case in Section 
4. (b) The ejected time of the 100 "Earths", r is the distance between the "Earths" and Kepler-16 A (and hereafter). The dotted line 
represents the location of Kepler-16 b. 
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Fig. 2. — "Earths" in Kepler-34, 35, 38 and 47 are represented in panels (a),(c),(e) and (g), respectively: The green filled circles represent 
stable habitable "Earths", while the red/blue circles stand for unstable/stable "Earths" outside the HZ after 10 Myr. The green and red 
dotted lines are the boundaries of the HZ and the initial locations of 100 " Earths" , respectively. The equilibrium temperature Tp and the 
location r of stable "Earths" in KepIer-34, 35, 38 and 47 are shown in panels (b),(d),(f) and (h), respectively. The error bars show the 
largest and smallest values of Tp and r. 
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Fig. 3. — Differences of (a) transit times (dTtr) and (b) TTV(dTTV) of detected planets with and without the perturbation of an 
additional "Earth". Because of the stronger perturbation between "Earth" and Kepler-16b in 1:1 orbital resonance, Ttr and dTTV of 
Keplcr-16b arc obviously larger than other four systems, and can be detected in three years. 



